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Non-monotonic magnetic field and density dependence of in-plane magnetoresistance 
in dilute two-dimensional holes in GaAs/AlGaAs 

H. Noh, Jongsoo Yoon,0 D. C. Tsui, and M. Shayegan 
Department of Electrical Engineering, Princeton University, Princeton, NJ 08544 

(February 1, 2008) 

We studied low temperature (T — 50mK) in-plane magnetoresistance of a dilute two-dimensional 
hole system in GaAs/AlGaAs heterostructure that exhibits an apparent metal-insulator transition. 
We found an anisotropic magnetoresistance, which changes dramatically at high in-plane fields 
(-B|l S^5T) as the hole density is varied. At high densities where the system behaves metallic at 
i3|| = 0, the transverse magnetoresistance is larger than the longitudinal magnetoresistance. With 
decreasing the hole density the difi'erence becomes progressively smaller, and at densities near the 
"critical" density and lower, the longitudinal magnetoresistance becomes larger than the transverse 
magnet oresist ance . 
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The scaling theory of localizationE predicts that, in the 
absence of electron-electron interaction, all states in two 
dimensions (2D) are localized in the zero temperature 
limit and that only an insulating behavior characterized 
by an increasing resistivity (p) with decreasing tempera- 
ture (T) is possible at low T. In contrast, experimental 
observations of a "metallic" behavior, characterized by 
a decreasing p with decreasing T, and thus an apparent 
metal-insulator transition (MIT) as charge carrier den- 
sity is lowered, are reported on many differeat low dis- 
order dilute 2D systems at low temperaturesJj To date, 
despite the large number of experimental as well as the- 
oretical papers in the literature, there is still no consen- 
sus on either the physics and the mechanism behind the 
metallic-like behavior or the origin of this apparent MIT. 

From the more recent experiments, two factors that 
strongly influence the charge transport close to this MIT 
in these low density high mobility 2D systems have 
emerged. First, disorder is found to play a significant 
role. A close examination of all such 2D systems re- 
veals that the critical density (ric for electrons and Pc 
for holes), the density above which the system shows the 
metallic behavior, is lower in a system with higher carrier 
mobility. The two most extensively studied, systems are 
the electron inversion layer in Si-MOSFETtljand the 2D 
holes in the GaAs/AlGaAs heterostructure.QB Typically, 
the 2D electrons in a high quality Si-MOSFET have a 
peak mobihty of ^ 2 x lO^cm^/Vs to 5 x lO'^cm^/Vs, 
and ric ^ 1 X 10^^ cm~^. In the extremely high mobility 
2D hole system (2DIIS) in GaAs/AlGaAs, which has a 
comparable effective mass at low densities and a more 
than 10 times higher peak mobility of ~ 7 x lO^cm^/Vs, 
Pc ^ Ix 10^° cm^^. Both Tic and Pc are found to decrease 
with decreasing disorder in the 2D system,Ll approaching 
the critical density characterized by the dimensionless 
density param|eter = 37, expected for Wigner crystal- 
lization in 2DJ3 

The second factor is the effect of a magnetic field ap- 
plied parallel to the 2D plane. The parallel field (i?||) 
is found to suppress the metallic behavior and also in- 
duce an unexpectedly large positive magnetoresistance 



(MR). In the case of the 2DHS in GaAsi this giant 
MR follows ~ exp(i?j|) at low and crosses into an 
~ exp{B^\) dependence at high B|| around a well defined 
characteristic field B* , which is hole density dependent. 
This behavior persists into the insulating phase, though 
B* becomes density independent. Similarly strong pos- 
itive MR was observed earlier in the electron inversion 
layer in Si-JVIOSFET's at low B|| though it is followed by 
saturations or in some cases a weak linear dependence at 
highBii. 

Recently, Das Sarma and HwangLl pointed out the im- 
portance of taking into account in the _B|| transport the 
finite thickness of the 2D electrons through which can 
be coupled into their orbital motion. They calculated 
the Boltzmann transport in the presence of iJy assuming 
(5-function random impurity scattering and a harmonic 
confinement of the 2D carriers to the heterojunction in- 
terface. They were able to produce a giant positive MR 
showing dependencies in qualitative agreement with 
experimental observations. In particular, their calcula- 
tion showed almost two orders of magnitude increase in 
resistance, when _B|| was increased to lOT. Furthermore, 
they found a crossover from a low field to high field be- 
havior due to the change of the relative magnitude of the 
cyclotron energy to the confinement energy with increas- 
ing field. The difference between the behavior seen in 
the Si inversion layer and that in the 2DIIS in GaAs in 
the high field regime was attributed to the different im- 
purity scattering mechanisms. Most importantly for ex- 
periments, they predicted that this giant in-plane MR is 
highly anisotropic: the transverse MR {p± , the MR when 
-B|| _L /, where I is the current) can be many times larger 
than the longitudinal MR (p||, the MR when iSy || /) for 
a moderately high in the lOT range. 

In the course of a systematic study on the nature of 
the insulating ground state in the extreme dilute limit of 
the 2DIIS in GaAs, we have investigated the By trans- 
port in one of the highest mobility low density sample 
currently available as a function of the hole density (p) 
across the apparent MIT with Bn up to 14T at dilution 
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refrigerator temperatures down to 50mK. By exploring a 
regime (lower density, lower temperature) not accessible 
in previous in-plane MR anisotropy experiments, we find 
a behavior which cannot be understood with currently 
available models. In addition to observing the expected 
behavior of a much larger MR in this regime, we find that 
the anisotropy, defined as the difference between p\\ and 
p± , changes sign in the high field region as the hole den- 
sity is decreased. For densities p > 2 x IO^"cni~^, is 
larger than p|| at high field similar to the model calcula- 
tion by Das Sarma and Hwang. However, the difference 
between and p\\ becomes progressively smaller as the 
density is decreased, and for densities p < 2 x 10^°cm~^, 
P|| becomes larger than pj^. In the rest of this paper, we 
describe in more detail the experiment and discuss our 
data in light of more recent results fropa. ■^He tempera- 
ture experiments in the metallic regime.E3 We also point 
out that the observed anisotropy is at variance with pre- 
dictions of existing theories explaining the large in-plane 
MR. 

Our samples are made from a modulation doped 
GaAs/AlGaAs (311)A single interface heterostructure 
grown by molecular beam epitaxy. The 2DHS at the 
(311)A interface even at B\\ — exhibits a mobil- 
ity anisotropjj known to be due to anisotropic surface 
morphologyO The mobility is high along the [233] crys- 
tallographic direction, and low along the [Oil] direction. 
We made two Hall bar (rectangular) samples, one for 
each direction, from the same wafer. For resistivity mea- 
surements, / was always passed along the sample length. 
Both samples also have back gates to change their den- 
sity. The peak mobilities at T = 50mK for our un- 
gated samples with a density of p = 4.2 x 10^°cm~^ are 
6.3_x lO^cmVVs along [233] and 4.2 x lO^cmVVs along 
[Oil]. These are among the highest mobilities achieved 
in samples of this type of single interface heterostruc- 
tures. Our measurements were made down to p = 
7.9 X lO^cm"^. This density is about 40% lower than that 
accessed by Papadakis et al. in their recent experiments 
using the GaAs/AlGaAs (311)A quantum wells at '^He 
temperatures.E3 It enabled us to measure the in-plane 
MR across the zero-field MIT at pc = 9.3 x lO^cm^^ p^j. 
a given density, we measured both p±{B\\) and p[|(i?||) 
by applying _B|| up to 14T. The samples were mounted, 
in separate cooldowns, on a rotating stage in a dilution 
refrigerator, which could adjust the angle between the 
2D plane and the direction of the magnetic field. They 
were mounted in such a way that either _L / for p± 
measurement or || / for p|| measurement, when the 
2D plane is aligned parallel to the magnetic field. Be- 
tween the pj_ and p|| measurements, the samples were 
thermally recycled, and measurements were repeated for 
consistency check. 

The in-plane MR, p±{B^\) and /9||(i?||), measured at 
50mK for / fixed in [233] direction arc shown in Fig. ^ (a) 
and (b), respectively, in the density range from 4.25 x 10^° 
cm~^ to 0.79 X 10^° cm~^. The sample shows insulating 
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FIG. 1. (a) p± and (b) py plotted as a function of -B|| at 
densities p = 4.25, 3.25, 2.86, 2.18, 1.56, 1.10, 1.00, 0.87, 0.79 
X 10^*^001"^ (from bottom to top). The critical density pc, 
where the MIT occurs at zero-field, is pc = 0.93 x 10^°cm~^. 
The characteristic field Be, above which the metallic behav- 
ior is suppressed, is marked by the short dashed line, and 
B* , which divides the low field and the high field regions, is 
marked by the long dashed line. The orientation of By and / 
are shown with respect to the crystal axes. In both (a) and 
(b), the solid lines are data measured by low frequency lock- in 
technique and the solid circles are those measured by DC I-V 
measurement, (c) Bc/B* plotted as a function of p. Solid 
(open) symbols are for By 1. I (By jj I). 



behavior at By = for the two lowest densities (the top 
two traces) and metallic behavior for the higher densities. 
Both px and py show strong increase with increasing By , 
about one order of magnitude for the highest density and 
more than two (p_l) or three (py) orders of magnitude for 
the lowest density at 14T. It is clear that the MR is larger 
for lower densities. We should also note that the MR in 
our sample is approximatelsi-ten times larger than that 
observed by Papadakis et aZ£3 at T = 0.3K in the density 
range where the two experiments overlap. 

For both pj^ and py , according to their dependence 
on _By , the p-B\\ plane can be divided into two regions: 
a low field region (By < B*) and a high field region 
(By > B*). In the low field region, p measured in both 
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field orientations increases rapidly with increasing 
following p ~ exp(i3y). This rapid increase slows down 
as B\\ is increased to the high field region, where p_L in- 
creases as exp[B\\), while p\\ shows a behavior depending 
on the density. For densities p > 2.18 x 10"'^°cm~^, also 
increases as exp{B\\), but for lower densities, it shows a 
stronger dependence, which becomes more pronounced 
when the density is lowered. For both B|| orientations, 
the metallic behavior observed at B\\ = and p > pc is 
suppressed for B\\ > Be (the short dashed lines in Fig. |l| 

(a) and (b)). Both Be and B* are smaller in Fig.Jl| (b) 
than in Fig. |l| (a). However, as shown in Fig. |^(c), 
when the ratio of Be to B* is plotted as a function of p, 
the two sets of data are nearly identical. These observa- 
tions result from a _B|| induced spin subbajid depopula- 
tion identified recently by Papadakis et a^JlJ and Tutuc 
et a/.E3 from their detailed analyses of the Shubnikov-de 
Haas oscillations observed at 0.3K. We shall later return 
to this point. 

In Fig. H (a), we plot both p^(_B|[) and from 
Fig. |l] normalized by the values at 5|| = in a single 
plot to show the anisotropy. Note that in Fig. (a) both 
sets of data were taken with / fixed along [233] and the 
direction of B\^ was varied from [Oil] to [233] to measure 
p± and p|| , respectively. In Fig. g (b), we show additional 
data, taken with / along [Oil] and -By along [233] 
and compare them with P|[_of Fig. |^ (a), which is taken 
with both / and By along [233]. In other words, in Fig. || 

(b) the direction of By is fixed along [233], and p± and 
/0|| data are shown for / along [Oil] and [233], respec- 
tively. The data of Fig. |2| highlight the central finding 
of our paper: in the high field region, on the right of the 
dotted lines, we observe an anisotropy in MR which re- 
verses sign as a function of density. At high densities 
p± is larger than py , while at low densities py is larger 
than p±. This py > p± anisotropy continues to grow 
across the MIT, and for p ^pc, in the B\\ — insulating 
side, p|| is approximately an order of magnitude larger 
than p± at 14T. Note that this trend is seen regardless 
of whether p^ and py are measured by fixing the direc- 
tion of / (Fig. I (a)) or B\\ (Fig. | (b)), suggesting that 
the relative orientation of / and _B|| is important in the 
high field region. No change of anisotropy at high field 
with decreaspg density was observed in earlier T = 0.3K 
experimentsjlj where a relatively high density sample, 
showing only the metallic behavior at B\\ — 0, was stud- 
ied. There, p_L was always larger than p|| even at -B|| up 
to 16T. Before we further discuss this surprising trend, 
it is useful to make a few remarks about the MR data at 
lower B\\ , on the left of the dotted lines in Fig. |2[ 

As mentioned before, each of our MR traces has a char- 
acteristic field, B* . It has been recently shown that B* 
marks tho-Seld above which the 2DHS becomes fully spin 
polarized Ji3'E3LJ Below B*, two spin subbands with un- 
equal populations are occupied while, above B*, there is 
only one occupied spin subband. The situation is simi- 
lar for 2D electron systems in Si-MOSFET's; there too a 
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FIG. 2. p± and py normalized by the zero-field resistivity 
p(0) and plotted together for (a) I fixed along [233] and (b) By 
fixed along [233]. The solid lines and solid circles are for p± 
and the dashed lines and open circles are for py . The densities 
for each pair of traces in (a) from bottom to top is p = 4.25, 
3.25, 2.86, 2.18, 1.56, 1.10, 1.00, 0.87, 0.79 xlO^°cm-^ The 
data for each p are offset by 10 as indicated by the scale bar 
for clarity. The densities in (b) from bottom to top is p = 
4.25, 2.18, 1.56, 1.10, 0.87 xlO^°cm"^ and the data are offset 
in such a way that zero-field values match with those in (a) 
for corresponding densities. The dotted lines are guide to the 
eye, representing the boundary between the low and high field 
regions. 

strong MR is observed at low fields when two spin sub- 
bands are populated, while the MR is essentially satu- 
rated once the system is fuilly spin polarised. oli3 Calcu- 
lations by Dolgopolov et alX^ and HerbutllZl attribute this 
low field MR to the change of screening as the spin polar- 
ization changes. In the case of 2DHS in GaAs (311)A, B* 
depends on the orientation of the i3y with respect to the 
crystal axes. In particular, B* is smaller for By 1|[233] 
compared to ]| [Oil] . This is due to the intrinsic energy 
band anisotropy of the GaAs (311)A 2DHS, as, iias-|bcen 
documented experimentally and thcorcticallyll3'll2Hl3 Our 
observation that B* in the pl data of Fig. ^ (a) is larger 
than in the py data of Fig. n (b) is consistent with these 
findings. That the data plotted as Bc/B* for the two -By 
orientations in Fig. ^ (c) should closely agree is also to be 
expected if the minority spin population is constant at 
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Be as found by Tutuc et alE3 and the spin polarization 
is linear in . 

It is likely that the MR anisotropy observed at low 
fields in Fig. H (a) is a consequence of the difference in 
spin subband depopulation rate for different orientations 
of with respect to the crystal axes. An alternative 
interpretation of the low field anisotropy may be based 
on a recent theory by Chen et uLIlB These authors find 
that an interplay between the spin-orbit coupling and 
the Zeeman splitting can lead to an anisotropic MR de- 
pending on the relative direction between / and . The 
anisotropy is small at low and high fields, and exhibits 
a maximum at an intermediate field when the Zeeman 
energy equals the spin-orbit induced spin-splitting. At 
first sight, this seems qualitatively consistent with the 
low field portion of Fig. |2| (a) where the difference be- 
tween p|| and p± goes through a maximum. However, 
this trend is not observed in Fig. |^ (b) which also com- 
pares and p± but now for a fixed direction. We 
conclude that the explanation of the anisotropic MR ob- 
served in the low field region in Fig. || (a) mostlikely lies 
in the anisotropic spin subband depopulationJlj 

We now return to a discussion of the high field data 
which are the main subject of this paper. We first note 
that these data are in a field region where only one spin 
subband is occupied and the system is fully spin polar- 
ized. In this region, calculations such as those reported 
in Refs. 16 and 17 predict that the MR would saturate. 
A saturation of the MR above B* is indeedLoften seen 
in 2D electron systems in Si-MOSFET'sEj-El This is in 
contrast, however, with the data shown in Fig. ^ where a 
large MR is observed even at high fields beyond B* . Cal- 
culations by Das Sarma and Hwangjj on the other hand, 
point to the coupling of to the orbital motion, and 
show that there can be a significant MR when the mag- 
netic length becomes comparable to or smaller than the 
layer thickness. They also predict that the MR should be 
anisotropic with p± > p|| . While this trend is seen in our 
data at high densities (Fig. §), the anisotropy is reversed 
at low densities. Their calculations predict that with 
decreasing density, the anisotropy {p± > p||) increases 
because of the increasing layer thickness. This is also 
at odds with our experimental data, although we should 
note that in our samples we expect the layer thickness 
to decrease with decreasing density since we are using a 
back gate. It is remarkable that the low density data in- 
dicate that at high fields p\\ > p± regardless of whether 
the direction of / or _B|| is kept fixed. This observation 
suggests that the anomalous anisotropy we observe at low 
densities may be intrinsic to very dilute 2D hole systems 
in GaAs. j—. 

Finally, Meirt£l explained the large in-plane MR in 
GaAs as from inhomogeneous distribution of charge car- 
riers due to disorder. In his model, the disorder poten- 
tial induces puddles of 2D holes which are interconnected 
by quantum point contacts (QPC), and the transport is 
dominated by the conduction through these QPC's. An 
in-plane magnetic field increases the bottom of the 2D 



subbands and effectively increases the conduction thresh- 
old energy of the QPC's, giving rise to the large positive 
MR. The MIT is a percolation transition which can be 
tuned by density and by in-plane magnetic field because 
the conduction through the QPC's depends on the rela- 
tive value of the Fermi energy to the threshold energy of 
the QPC's. However, no anisotropy is expected in this 
model, where the MR arises from the shift of the 2D 
subband energy by the presence of . 

In summary, we have measured the in-plane MR p± 
and of a very high mobility, dilute 2DHS in the 
GaAs/AlGaAs heterostructure in the density range from 
4.2 X lO^^cm^^ to 7.9 x lO^cm"^, spanning across the ap- 
parent zero-field MIT. In the low _B|| region, the observed 
anisotropy is consistent with the spin subband depopu- 
lation which depends on the orientation of -B|| relative to 
the crystal axes. In the high field region, we observed a 
change in the sign of the anisotropy with decreasing den- 
sity. At high densities, p± > p\\. However, as the den- 
sity is lowered, this anisotropy decreases and p|| becomes 
equal to p± around ~ 2 x 10"'^°cm~^. For still lower den- 
sities, p|| > p±, and the anisotropy changes sign and also 
grows continuously across the MIT. These experimental 
observations are not explicable with existing theoretical 
models for the in-plane MR. 
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